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Abstract 

A popular aspect of black holes physics is the mathematical analogy between their laws, coming 
from general relativity and the laws of thermodynamics. The analogy is achieved by identifying 
a suitable set of observables, precisely: (a) E = M (being E the thermodynamic free energy and 
M the mass of the BH), (h) T = ok (with T the absolute temperature, k the so-called surface 
gravity on event horizon and a a suitable dimensional constant) and (c) S = (l/87rQ!)yl (where S 
is the thermodynamic entropy of the black hole and A the surface of the event horizon). However, 
despite numerous investigations and efforts spent on the subject, the theoretical foundations of 
such identifications between physical quantities belonging to apparently unrelated frameworks are 
not yet clear. The goal of this work is to provide the contribution to the black hole entropy, coming 
from matter in the black hole exterior. We propose a classical solution for the kinetic description 
of matter falling into a black hole, which permits to evaluate both the kinetic entropy and the 
entropy production rate of classical infalling matter at the event horizon. The formulation is based 
on a relativistic kinetic description for classical particles in the presence of an event horizon. An 
H-theorem is established which holds for arbitrary models of black holes and is valid also in the 
presence of contracting event horizons. 

PACS numbers: 51.50+v, 52.20-j, 52.27.Gr 
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I. INTRODUCTION 



The remarkable mathematical analogy between the laws of thermodynamics and black 
hole (BH) physics following from classical general relativity still escapes a complete and 
satisfactory interpretation. In particular it is not yet clear whether this analogy is merely 
formal or leads to an actual identification of physical quantities belonging to apparently 
unrelated framework. The analogous quantities are E <-> M, T ^ an and S ^ (l/87ra;)y4, 
where A and k are the area and the surface gravity of the BH, while a is a constant. 
A immediate hint to believe in the thermodynamical nature of BH comes from the first 
analogy which actually regards a unique physical quantity: the total energy. However, at 
the classical level there are obstacles to interpret the surface gravity as the BH temperature 
since a perfectly absorbing medium, discrete or continuum, which is by definition unable 
to emit anything, cannot have a temperature different from absolute zero. A reconciliation 
was partially achieved by in 1975 by Hawking [l], who showed, in terms of quantum particle 
pairs nucleation, the existence of a thermal flux of radiation emitted from the BH with a 
black body spectrum at temperature T = hK,/27TkB {Hawking BH radiation model). The 
last analogy results the most intriguing, since the area A should essentially be the logarithm 
of the number of microscopic states compatible with the observed macroscopic state of the 
BH, if we identify it with the Boltzmann definition. In such a context, a complication arises 
when one strictly refers to the internal microstates of the BH, since for the infinite red shift 
they are inaccessible to an external observer. An additional difficulty with the identification 
S ^ {l/8TTa)A, however, follows from the BH radiation model, since it predicts the existence 
of contracting BH for which the radius of the BH may actually decrease. To resolve this 
difficulty a modified constitutive equation for the entropy was postulated , in order to 
include the contribution of the matter in the BH exterior, by setting 

1 A 

S'-S^-k-, (1) 



(S" denoting the so-called Bekenstein entropy) where S is the entropy carried by the matter 

4'^ Gh 



outside the BH and Sbh = identifies the contribution of the BH. As a consequence a 



generalized second law 

6S' > (2) 
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was proposed which can be viewed as nothing more than the ordinary second law of 

thermodynamics apphed to a system containing a BH. From this point of view one notices 
that, by assumption and in contrast to the first term S, Sth cannot be interpreted, in a 
proper sense, as a physical entropy of the BH, since, as indicated above, it may decrease in 
time. 

This approach however is unconvincing since the precise definition and underlying sta- 
tistical basis both for S and Sbh remain obscure. Thus a fundamental problem still appears 
their precise estimates based on suitable microscopic models. Since the evaluation of Sbh 
requires the knowledge of the internal structure of the event horizon (excluding for causality 
the BH interior), the issue can be resolved only in the context of a consistent formulation of 
quantum theory of gravitation [4, |5|. This can be based, for example, on string theory [6|] 
and can be conveniently tested in the framework of semiclassical gravity [?, 8|. Regarding, 
instead the entropy produced by external matter S, its evaluation depends on the nature, 
not yet known, of the BH. However, even if one regards the BH as a purely classical object 
surrounded by a suitably large number of classical particles its estimate should be achievable 
in the context of classical statistical mechanics. 

In statistical mechanics the "disorder" characterizing a physical system, classical or quan- 
tal, endowed by a large number of permissible microstates, is sometimes conventionally 
measured in terms of the so-called Boltzmann entropy 5*^ = Kin W. Here K is the Boltz- 
mann constant while is a suitable real number to be identified with the total number 
of microscopic complexions compatible with the macroscopic state of the system, a number 
which generally depends on the specific micromodel of the system. Therefore, paradoxi- 
cally, the concept of Boltzmann entropy does not rely on a true statistical description of 
physical systems, but only on the classification of the internal microstates (quantal or clas- 
sical). As is well known[10], Sb can be axiomatically defined, demanding (i) that it results 
a monotonic increasing function of W and (ii) that it satisfies the entropy additivity law 
Sb{WiW2) = 5*^(^1) + Sb(W2). Boltzmann entropy plays a crucial role in thermodynamics 
where (i) and (ii) have their corresponding laws in the entropy nondecreasing monotonicity 
and additivity. Since in statistical mechanics of finite system it is impossible to satisfy both 
laws exactly, the definition of Sb is actually conditioned by the requirement of considering 
systems with W >> 1 (large physical systems). 

An alternate definition of entropy in statistical mechanics is the one given by the Gibbs 
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entropy, in turn related to the concept of Shannon information entropy. In contrast to the 
Boltzmann entropy, this is based on a statistical description of physical systems and is defined 
in terms of the probability distribution of the observable microstates of the system. In many 
cases it is sufficient for this purpose to formulate a kinetic description, and a corresponding 
kinetic entropy, both based on the one-particle kinetic distribution function. In particular, 
this is the case of classical many-particle systems, consisting of weakly interacting ultra 
relativistic point particles, such as those which may characterize the distribution of matter 
in the immediate vicinity of the BH exterior. 

These issue have motivated a recent research effort on the subject by the present authors 
91 . The primary goal has been, in particular, to provide an explicit expression for the 
contribution S, which characterizes Bekenstein law ([T]), to be evaluated in terms of a suitable 
kinetic entropy, and to estimate the corresponding entropy production rate due to infalling 
matter at the BH event horizon. In addition we intend to establish an H-theorem for the 
kinetic entropy which holds, in principle, for a classical BH characterized by event horizons 
of arbitrary shape and size and even in the presence of BH implosions or slow contractions. 
This is obtained in the framework where the classical description of outside matter and space 
is a good approximation to the underlying physics. 



II. KINETIC DESCRIPTION OF INFALLING MATTER 

The basic assumption is that the matter falling into the BH is formed by a system Sjsi 
of classical point particles moving in a classical spacetime and described by Hamiltonian 
dynamics, while the event horizon can be treated as a classical absorbing porous wall. We 
adopt for this purpose a covariant kinetic formalism taking into account the presence of an 
event horizon. The evolution of such a system is well known and results uniquely determined 
by the classical equations of motion, defined with respect to an arbitrary observer O. To 
this purpose let us choose O, without loss of generality, in a region where space time is 
(asymptotically) flat, endowing with the proper time r, with r assumed to span the set / C i? 
(observer's time axis). Each particle is described by the canonical state x, which in the case 
of point particles spans the 8— dimensional phase space F, where x = {r'^,Pn). Moreover, 
its evolution is prescribed in terms of a suitable relativistic Hamiltonian H = H (x) so that 
the canonical state x = {r^,Pn) results parameterized in terms of the world line arc length 
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s (see [Ul]). As a consequence, requiring that s = s{t) results a strictly monotonic function 
it follows that, the particle state can be also parameterized in terms of the observer's time 
r. To obtain the a kinetic description for such a system we require, as usual, ^ 1. In 
addition, it is assumed that interactions between point particles of take place only via a 
mean-field Hamiltonian H (x) and hence that S^- can be identified with a weakly interacting 
relativistic gas. For Sn we introduce the kinetic distribution function for the observer O, 
defined as follows 



Pg(x) = p(x)5(s - s(r))5( - 1) (3) 

where p (x) is the conventional kinetic distribution function in the 8— dimensional phase 
space, to be assumed suitably smooth and summable in F {Assumption a of regularity). 
Notice that the Dirac deltas here introduced must be intended as physical realizability equa- 
tions. In particular the condition placed on the arc length s implies that the particle of the 
system is parameterized with respect to s(r), i.e., it results functionally dependent on the 
proper time of the observer; instead the constraints placed on 4-velocity implies that is a 
tangent vector to a timelike geodesic. The event horizon of a classical BH is defined by the 
surface vh specified by the equation 

R{x) = th (4) 

where x denotes a point of the space time manifold, while R{x) reduces to the radial coor- 
dinate in the spherically symmetric case. 

According to a classical point of view, let us now assume that the particles are "captured" 
by the BH (i.e., for example, they effectively disappear for the observer since their signals 
are red shifted in such a way that they cannot be anymore detected [l^) when they reach 
a suitable surface 7 {capture surface) of equation 

Re{x) = r,. (5) 

Here = (1 + e)rj^, while e > depends on the detector and the 4— momentum of the 
particle. The presence of the BH event horizon is taken into account by defining suitable 
boundary conditions for the kinetic distribution function on the hypersurface 7, to be treated 
as an effective absorbing porous wall. For this purpose, we distinguish between incoming 
and outgoing distributions on 7, Pg(x) and Pg(x) corresponding respectively to UaU" < 



and HaU" > 0, where locally radial outward 4— vector. Therefore, the boundary 

conditions on 7 are specified as follows 

Pg W = pW^(^ - sirrn,/^ - 1) (6) 
p+(x) = 0. (7) 

It follows that it is possible to represent the kinetic distribution function in the whole space 
time manifold in the form 

Pg (x) = Pg (x) + p^ (x) (8) 

where 

p±(x) = p(x)5(s-s(r))5(^/^-l) X 

xe±(i?e(x)-r,(s(r))) (9) 

with respectively denoting the strong and the weak Heaviside functions 

i 1 for a>Q 
e-{a) = { - (10) 

I for a < 0. 

and 

^ _ I 1 for a> ^ ^ 

0+(a) = (11) 

I /or a < 0. 

In the sequel we shall introduce the hypothesis that the distribution P(5(x) has a compact 
support {Assumption [3). 

It is important to stress that in the definition of the boundary conditions no detailed 
physical model is actually introduced for the particle loss mechanism, since all particles are 
assumed to be captured on the same hypersurface 7, independent of their mass, charge and 
state. This provides a classical loss model for the BH event horizon. 

Let us now consider the evolution of the kinetic distribution function Pg(x) in external 
domains, i.e. outside the event horizon. Assuming that binary collisions are negligible, 
or can be described by means of a mean field, and provided that the phase space volume 
element is conserved, it follows the collisionless Boltzmann equation, or the Vlasov equation 
in the case of charged particles j2o|, 

ds fdrf'dpix) ^ dp^dp{x) \ _^ ^^^^ 



dr { ds dri^ ds dp^ 



with summation understood over repeated indexes, while p(x) denotes evaluated at 



dro{s) 



. This equation resumes the conservation of the 

s=s(t) 



= r°(s(r)) and po = m 
probability in the relativistic phase space in the domain external to the event horizon. 
Invoking the Hamiltonian dynamics for the system of particles, the kinetic equation takes 
the conservative form 

-[p(x),i^L = 0. (13) 

III. BH STATISTICAL ENTROPY AND H-THEOREM 

Let us now introduce the appropriate definition of kinetic entropy S{p) in the context of 
relativistic kinetic theory. We intend to prove that in the presence of the BH event horizon 
it satisfies an H theorem. 



The concept of entropy in relativistic kinetic theory can be 
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brmulated by direct extension 



15|, 



16| . For this purpose we 



of customary definition given in nonrelativistic setting 
introduce the notion of kinetic entropy, measured with respect to an observer endowed with 
proper time r, as follows 

S{p) = -pj d^{s)6{s - sirM^/^- l)p(x) Inp(x), (14) 
r 

where p(x) is strictly positive and, in the 8-dimensional integral, the state vector x is param- 
eterized with respect to s, with s denoting an arbitrary arc length. Here P is the principal 
value of the integral introduce in order to exclude from the integration domain the subset 
in which the distribution function vanishes. Hence S{p) can also be written: 

S{p) = -pj d^{s)6{s - s(r))pi(x) Inp(x), (15) 
r 

where pi(x) now reads 



pi(x) = e(r(.) - r,(.))5( - l)p(x(s)). (16) 

In the sequel S{p) will be denoted as BH classical entropy. Differentiating with respect to 
T and introducing the invariant volume element d^rd^p, the entropy production rate results 
manifestly proportional to the area A of the event horizon and reads 

^-f^ + ^ = -p/ dhd'pF,.^ [6 (r - r,) plnp] , (17) 



where Frr^ is the characteristic integrating factor 



_ dsjr) ( dr dr,. 



dr \ds ds 

Indeed, the r.h.s represents the entropy flux across the event horizon while ^ and ^ are 
the contributions to the entropy production rate which depend, respectively, on the velocity 
of the infalling matter and of the event horizon: 

"^^^ - P [dhd'p^^6{r-r,)p\np, (19) 



dr J dr ds 

r- 

- P [dhd'p^^^5ir-r,)p\np. (20) 



dr J dr ds 

r- 

As a consequence, ^ and ^ can be interpreted as the contributions to the BH entropy 
production rate carried respectively by the infalling matter and the event horizon. Here, F" 
is the sub-domain of phase space corresponding to the particle falling into the BH. Hence, 
it follows that in the above integrals,^ results by definition negative. Instead, ^ has not 
a definite sign and it can be negative in the case of contracting event horizons. We can 
also write the above expression in terms of the kinetic probability density evaluated at the 



hypersurface ^Uy,u^ — 1, defined as /(x) = pjN. It follows plnp = NflnNf. At this 
point we adopt a customary procedure in statistical mechanics 17|] invoking the inequality 



NflnNf > Nf -1 (21) 

and notice that in the sub-domain be in which Frr^ > there results by definition p = 0. 
Hence it follows that in F" 

Frr, < (22) 

where by construction > 0. This result holds independent of the value of Next, 
due to Assumption a and (3 on the surface 7 the kinetic distribution function results nec- 
essarily non-zero only in a bounded subset of phase space. Therefore, if 6 is an arbitrary 
infinitesimal, there exists a bounded subset C F~ such that Nf results infinitesimal (of 
order 6) in the complementary set F~ \ i7 (which includes the set of improper points of F~) 
and moreover 

P j dhd^p\Frr^S{r -r,) NflnNf ~ 0(5 In (5). (23) 

r-\n 



8 



At least in a subset of Q, Nf is by assumption positive and such that Nf > S. Therefore 
there results 

0<p/A.3p|F„.|M.-r.).M (24) 

n 

where M is a suitable finite constant. Thus one obtains 
dSip) 



dr 



>pj dhd^p \Frr, I 5 (r - rj [Nf - 1] + O (5 In 5) (25) 



n 



The first term of the r.h.s of fl25|) can be interpreted in terms of the effective radial velocity 
of incoming particles 

V;ff ^-[d'p \F„, \5ir- r,) /, (26) 



no, 
( 

while no is the surface number density of the incoming particle distribution function. Finally 
we invoke the majorization 

^ >^>NMl[dhnoV;ff' 
dr dr [J^ 

-M + 0{6\n6) (27) 



and notice that, due to the arbitrariness of 6, N can always be chosen sufficiently sufficiently 
large to satisfy the inequality 5 > 0. We stress that inf {J d^vnoV^^^^ can be assumed 
strictly positive for non isolated BH's surrounded by matter while 6 is arbitrary. This yields 
the relativistic H-theorem: 

THEOREM - H- Theorem for the BH classical entropy 

In validity of the kinetic equation l[T^) . boundary conditions ^ and regularity con- 
ditions defined by Assumptions a and f3, the BH classical entropy S{p) results uniquely 
defined, according to Eq.^^. Moreover, if the number N ^ 1 of classical point particles of 
the system Sn is assumed suitably large, there results: 

dS{p) ^dSi^dS2 ^ Q ^28) 
dr dr dr 2 ~ ' 
where ^ and ^ are respectively the contributions to entropy production rate ( I7P|) . due 

to respectively to infalling matter and the event horizon. In particular there results ^^^^ = 0, 

if and only if uq = 0, being uq the number density on the boundary surface 7. Moreover, 

the inequality holds even if: 



dVf , , 



(contracting event horizon). 
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IV. CONCLUSIONS 



Let us briefly analyze the basic implications of this result. First we notice that the H- 
theorem here obtained appears of general validity even if achieved in the classical framework 
and under the customary requirement ^ 1 (large classical system). Indeed the result 
applies to BH having, in principle, arbitrary shape of the event horizon. The des crip tion 
adopted is purely classical both for the falling particles (charged or neutral 18|, ll9|, |20|, l2l| ) 
and for the gravitational fleld and is based on the relativistic coUisionless Boltzmann equation 
and/or the Vlasov equation respectively for neutral and charged particles. A key aspect of 
our formalism is the deflnition of suitable boundary conditions for the kinetic distribution 
function in order to take into account the presence of the event horizon. The expressions 
for the entropy and entropy production rate, respectively Eqs. ffT^ and f|T7|) . can be used for 
speciflc applications to astrophysical problems. Finally interesting features of the derivation 
are that the entropy production rate results proportional to the area of the event horizon and 
that the formalism is independent of the detailed model adopted for the BH. In particular 
also the possible presence of an imploding star (contracting event horizon) is permitted. 
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